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ABSTRACT 

This paper uses numerical simulations to test the formation time distribution of dark 
matter haloes predicted by the analytic excursion set approaches. The formation time 
distribution is closely linked to the conditional mass function and this test is there- 
fore an indirect probe of this distribution. The excursion set models tested are the 
extended Press-Schechter (EPS) model, the ellipsoidal collapse (EC) model, and the 
non-spherical collapse boundary (NCB) model. Three sets of simulations (6 realiza- 
tions) have been used to investigate the halo formation time distribution for halo 
masses ranging from dwarf-galaxy like haloes (M = 10~ 3 M*, where M* is the charac- 
teristic non-linear mass scale) to massive haloes of M — 8.7M*. None of the models 
can match the simulation results at both high and low redshift. In particular, dark 
matter haloes formed generally earlier in our simulations than predicted by the EPS 
model. This discrepancy might help explain why semi-analytic models of galaxy for- 
mation, based on EPS merger trees, under-predict the number of high redshift galaxies 
compared with recent observations. 

Key words: methods: N-body simulations - cosmology:theory -dark matter - 
galaxies:haloes - galaxies: formation 



1 INTRODUCTION 

In the present preferred cosmological model, dark matter 
haloes form hierarchically through accretion and merging of 
smaller structures that grow from a Gaussian initial den- 
sity field. This process is modeled by Press-Schechter the- 
ory! Press & Schechter 19741, which simply assumes a re- 
gion with the mass over-density above a certain threshold 
will turn around and eventually collapse to form a bound 
object. Bond et al.(1991) developed an excursion set ap- 
proach and used it to derive the number density of collapsed 
dark matter haloes more rigorously. Lacey & Cole (1993) 
used the excursion set theory to predict the merger rate at 
which small objects merge with each other to form larger ob- 
jects, the conditional probability of its progenitor for a par- 
ent halo, as well as the survival probability of haloes. The 
predictions for the halo formation time distribution have 
been tested against N-body simulations [with 128 3 particles] 
l |Lacey fe Cole 1 994 1 . The conditional probability predicted 
by the excursion set theory (or EPS theory) has been widely 
used to plant merger trees (jKauffmann & White ~993l 

* E-mail address: linwp@center.shao.ac.cn, 

ypjing@center.shao.ac.cn, d90222005@ms90.ntu.edu.tw 



|Lacey fc Cole 1994|ISomerville fc Kolatt 19991 so as to con- 
struct semi-analytical models of galaxy formation, and to 
study the clustering of dark matter halos (Mo & White 1996; 
Mo, Jing, & White 1997). However the EPS theory may fail 
to model the details of halo formation. For instance, the 
predicted mass function was not found to fit the simulation 
results very well (e.g. Lee & Shandarin 1998, Sheth & Tor- 
men 1999). Tormen (1998) also reported that the excursion 
set predictions did not well fit the conditional mass function 
of sub-clumps in simulations. In addition, since the EPS 
theory failed to correctly describe the spatial distribution of 
small haloes in high resolution numerical simulations (Jing 
1998, 1999; Lee & Shandarin 1999; Sheth & Tormen 1999), it 
may also fail to predict the halo formation time distribution. 
There are already some indications that the distribution of 
halo formation time in N-body simulations is not consistent 
with predictions llvan den Bosch 2002a). 

In fact, it is an approximation that the spheri- 
cal collapse assumed in the EPS formalism. This pro- 
cess can be better modeled by a triaxial turn-around 
model (|Bond fc Myers 1996| |Sheth, Mo fc Tormen 200T| 
ISneth fe~^brrrien^00"2T ~ One of such models is the so-called 
Ellipsoidal Collapse model or EC model. By taking the ellip- 
soidal collapse into account, Sheth, Mo and Tormen (2001) 
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found that the modified mass function can fit N-body sim- 
ulations well. Sheth & Tormen (2002) pointed out that the 
conditional mass function was not universal since it was not 
consistent with their simulation results at every redshift. An- 
other model is the Non-spherical Collapse Boundary model 
(or NCB model shortly) proposed by Chiueh & Lee (2001). 
It relates the halo formation to the collapse of the Zel'dovich 
pancakes. A recent version of the model was presented by 
Lin, Chiueh & Lee (2002). Both the EC and the NCB mod- 
els were calibrated to fit the spatial two-point correlation 
function of halos (Jing 1998) and the mass function (Sheth 
& Tormen 1999) over a large range of halo mass. 

This paper uses three sets of high resolution N-body 
simulations [with 256 3 or 512 3 particles] to study the distri- 
bution of halo formation redshift and make comparison with 
theory predictions. It is organized as follows. In section |5] 
we present the theory predictions by three analytical models. 
The simulations are described briefly in section together 
with the method to find the halo formation redshift. There 
the simulation results are compared to the predictions. Main 
conclusions and discussion are given in section 0] 



2 THE THEORETICAL PREDICTIONS FOR 
THE DISTRIBUTION OF HALO 
FORMATION REDSHIFT 

The halo formation redshift is defined as the redshift at 
which its main progenitor has accumulated half of the halo 
mass. According to the EPS theory, the probability that a 
volume of mass Mi , which is within the region of mass M2 
collapsed at redshift Z2, collapsed to form a progenitor at 
redshift zi is given by the conditional probability function, 
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where S c i = S c (zi) and <5 C 2 = 5 c (zi) l |Lacey fc Cole 1993| . 
Here S c (z) is the critical over-density required for the spher- 
ical collapse at redshift z (see eq.A6), and Si = a 2 (Mi) 
where a (Mi) is the rms of the initial density fluctuation field 
smoothed on a scale which contains mass Mi, extrapolated 
using linear theory to the present time (see eq. Al). 

Thus the conditional probability can be converted into 
the probability that a halo of mass M2 at redshift Z2 has a 
progenitor in the mass range Mi to Mi + dMi at an earlier 
epoch z\ (i.e. the conditional mass function), 

dp 
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When using the excursion set model with ellipsoidal 
collapse (Sheth & Tormen 2002J, the conditional probabil- 
ity should be replaced by following function, 
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where the moving barrier B(S,z) = VaS» [1 + P(S/aS,) a ] 
with S* = S 2 (z). The parameters are adopted from the best 
fitting of the mass function with N-body simulations, a = 
0.707, a = 0.485, j3 = 0.615 jSheth, Mo fc Tormen 2001] 

ISheth fc Tormen 20021 . 

While for the NCB model HChiueh fc Lee 20011, the 
conditional probability is a fitting formula IILin et al. 2002 1 
which reads as, 
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with x = (\/S2/5 c (z2) — 0.25)/k, and q(n) is required to 
satisfy the normalization condition 
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Here k is the separation of two boundaries, which is defined 
as S c (zi)/S c (z 2 )- Refer to the appendix of Lin et al. (2002) 
for further explanation of the fitting procedure for these pa- 
rameters. Note that we use the uncorrected [eq.(29) of 
Lin et al.] rather than the corrected jj, [eq.(33) of Lin et 
al.]. There are two reasons for this choice. First, the halo 
formation distribution is incorrectly normalized when using 
the corrected [i . Second, the correction for // could be mass 
dependent. 

Integrating equation (2) over the mass range M2/2 < 
Mi < M2 gives the probability P(< ti) that its formation 
time is earlier than ti or its formation redshift is larger than 
z\. For a halo with mass Mo formed at the present, we set 
M2 = Mo,t2 = to,Z2 — zo = 0. Therefore the probability 
can be written as 



P(< t f ) = P(> z f ) 



So 
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with S h = S(M Q /2) | |Lacey fc Cole 1993> . In this integra- 
tion, the conditional probability functions (1), (4) and (5) 
are used for the EPS, EC and NCB models respectively. The 
accumulative probability P(> Zf) and thereby its redshift 
distribution can be calculated numerically. 

The predictions of the halo formation redshift distribu- 
tion for 3 typical masses are shown in Fig0 Here we assume 
a LCDM cosmogony (with as = 0.9) whose parameters will 
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Figure 1. The model predictions of the formation rcdshift for 
haloes with mass of 10 — 3 , 0.1 and 10M*. The left panels give the 
differential probability distribution of the halo formation redshift, 
and the right panels give the accumulative probability of haloes 
formed at redshift higher than z<. 



be given in next section. The left and right panels show re- 
spectively the formation redshift distribution and the prob- 
ability that a halo formed at redshift larger than Zf. The 
solid, dotted and dashed lines represent the EPS, EC and 
NCB predictions respectively. As can be seen, for haloes with 
a small mass of 10 _3 M, , the predictions of the three mod- 
els differ dramatically from each other. At low redshift the 
NCB prediction is close to the EPS one, however the pre- 
dictions by the NCB and EC model are coincident at high 
redshift. For haloes with masses of 0.1M, and 10M,, the 
EC and NCB results only have a small difference. For the 
large mass, there is almost no big difference in the prediction 
among the three models. In general, the dP/dzj profiles of 
the EC and NCB models are broader and have lower peaks 
than those of the EPS model. Compared to the EPS model, 
the EC or the NCB model predicts a larger fraction of haloes 
formed at high redshift and the EC model predicts a smaller 
fraction at low redshift. 



3 SIMULATIONS 

Three samples of N-body simulations are used to study the 
formation redshifts of dark haloes with mass ranging from 
10 -3 M* to 8.7 M*. Each simulation has at least 30 outputs, 
and can be used to trace the formation of a halo accurately. 
The cosmological model is the currently popular flat low- 
density model with the density parameter flo = 0.3 and the 
cosmological constant Ao = 0.7 (LCDM). The shape pa- 
rameter of the linear density power spectrum is V = Qoh = 
0.2. The characteristic mass M, is 9.55 x 10 12 /i _1 M Q for 
LCDMa, LCDMc simulations and 1.66 x 10 13 /i _1 M P> for 



LCDMb simulations. Other parameters of the simulations 
are listed in Tabled where as is the amplitude of the power 
spectrum, TV is the total number of particles in the simula- 
tion boxes, m p is the mass of a particle and z% is the initial 
redshift of the simulations. 

The simulation data were generated on the VPP5000 
Fujitsu supercomputer of the national Astronomical Obser- 
vatory of Japan with a vectorized-parallel P 3 M code (Jing 
& Suto 1998, Jing & Suto 2002). 



3.1 The formation of small dark matter haloes 

We use the LCDMa simulations to study the formation of 
small dark haloes. For each simulation, 56 out of the 169 
outputs are used. Only those haloes with 100 particles or 
more at the present are included to assure a reliable identi- 
fication of haloes (i.e., at the half-mass formation redshift, 
they have 50 particles at least). The haloes with more than 
1200 particles will not be considered, since these haloes are 
not abundant enough to have a reliable determination of the 
distribution of the halo formation redshift. Thus the dark 
haloes we study here have a mass range between 7.73 x 10 9 
and 9.28 x 10 10 h' 1 M© (~ 10~ 3 - 10~ 2 M»). These haloes 
are "small" and represent the typical dark haloes of dwarf 
galaxies. 

Dark matter haloes are identified with the spherical 
over-density method ( |Jing fc S uto 2002} . We adopt the fol- 
lowing method to find the formation redshift of a halo. At 
the beginning we pick up a halo as a parent halo in the final 
output (2 = 0), and find out the particles within its virial 
radius. Then the member particles are traced in the outputs 
from high redshift to low redshift step by step. We calculate 
the fraction of the member particles in all progenitors, and 
select the progenitor with the maximum number of parti- 
cles as the main progenitor. Generally, the members of the 
main progenitor will increase with time because of merg- 
ing and accretion, although in a few circumstances its mass 
may decrease because of mass loss due to the tidal stripping 
by nearby haloes and/or unbound particles. When the main 
progenitor has reached at least half of the parent halo's mass 
for the first time, we define the corresponding redshift as the 
formation redshift (half-mass formation redshift zj) of the 
selected parent. Alternatively, we can go along the merger 
tree from redshift to high redshift, and we can define the 
formation redshift as the time when the mass of the main 
progenitor first drops below half of the parent halo's mass. 
We find that these two definitions give almost an identi- 
cal distribution of the halo formation time distribution 1 , as 
serious tidal stripping of the main progenitor happens only 
rarely. In the following, we will adopt the first definition and 
compare the simulations with the three analytical models. 

The results are shown in Fig. IHa)-(c). The solid points 
are the result measured from the simulations, and the verti- 
cal error bars present the la error of the mean value derived 
from the scatter between the realizations. The solid, short- 
dashed, and long-dashed lines are the predictions of the EPS, 
EC and NCB analytical models. For the model predictions, 
we plotted two lines for halo at the lower mass end (the 



1 There is only a slightly shift of the profile, but the conclusions 
change little. 
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Table 1. Model parameters for simulations 



Model 
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h- 1 Mq 
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7.73 x 10 7 
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1200 


36 


60 


1 
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Figure 2. The distribution of the halo formation redshift for 
dwarf-galaxy like haloes. The points with error bars are the sim- 
ulation results, and the lines are the predictions of the three ana- 
lytical models. Different panels show the results for haloes of N p 
particles. The total number of simulated haloes (2 realizations) 
used for the analysis are 2334/2295, 1747/1786, 505/456 for the 
results in panel (a)-(c) respectively. 



redshift should fall between this output and the earlier one. 
This uncertainty is shown by the right-hand "error" bars. 

As the figure shows, at low redshift (< 1) the EPS and 
NCB models predict more halos than the simulations, while 
the EC predictions are consistent with the simulation re- 
sults. However at high redshift, the EPS model fits the sim- 
ulation results better than the EC and NCB models. In every 
panel, the formation redshift distribution of the simulations 
peaks at a higher redshift than the EPS prediction 2 and 
has a narrower profile (but with a higher peak) than the EC 
and NCB predictions. 

3.2 The formation of sub-A/„ dark matter haloes 

The LCDMb simulations are used to study the formation 
of sub-Af* dark haloes. Again the number of particles in a 
halo spans from 100 to 1200, or the corresponding halo mass 
ranges from 4.95 x 10 11 h' 1 M e to 5.94 x 10 12 h' 1 M (i.e. 
~ 0.03— 0.36M* ). The definition of haloes is slightly different 
from that used in the last subsection. Here distinct haloes 
were found using the FOF method with a bonding length 
0.2 times of the mean particle separation. We have tested 
the results for one realization against the two halo identi- 
fication algorithms, and found that the two identifications 
give nearly identical results. We calculate the halo formation 
redshift as in the last subsection for the 30 outputs, and plot 
the results in Fig0 

The results found for the sub-M* halos are quite sim- 
ilar to those for smaller haloes, but we can see continuous 
changes with halo mass of the model predictions relative to 
the simulation results. The EPS model predicts too many 
haloes of low formation redshifts again, while the EC and 
NCB models predict relatively well the fraction of these 
haloes. On the other hand, the EC and NCB models pre- 
dict too many haloes of high formation redshift. Note that 
for the mass range chosen here the EC and NCB predictions 
are close to each other (cf. Fig0. 



curve with higher peak) and the higher mass end (the curve 
with lower peak) respectively. We plot both of them, be- 
cause there may be a bias when using the prediction for the 
mean halo mass, as the formation probability is invariant in 
terms of a (M) 2 (rather than M). For the narrow ranges of 
halo mass chosen in this study, this bias is sufficiently small 
as the predictions at the upper and lower mass limits are 
very close. To elucidate another possible bias caused by the 
limited output resolution, we also add a right-hand "error" 
bar for the simulation results. This bias is introduced by the 
fact the halo formation time in simulation is measured to be 
the redshift at which the fraction of the member particles is 
larger than 0.5 for the first time, however the real formation 



3.3 The formation of large dark matter haloes 

We use a simulation (LCDMc) with 512 3 particles (Jing 
2002) to study the formation of large dark haloes. The cor- 
responding halo mass ranges from 1.67 x 10 12 h' 1 M to 
8.35 x 10 13 h' 1 Mq (~ 0.17 - 8.74M*). We did the same 

2 For the theoretical predictions of the halo formation redshift, 
a test on the infrared cut-off of wavenumber k due to a finite 
simulation box has been done. The cut-off can modify the M — 
<t 2 (M) relation so that it can change the distribution of the halo 
formation redshift. However, for the small haloes considered here, 
the wavenumber cutoff has negligible effect on their formation 
redshift. 
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Figure 3. The same as Fig[5] but for sub-Ai* haloes. The total 
number of simulated haloes (3 realizations) used for the analysis 
are 3577/3487/3607, 2698/2549/2619, 758/773/724 for the results 
in panel (a)-(c) respectively. 



Figure 4. The same as Fig|5| but with 512 3 particles in the 
simulation box and for haloes with mass > 0.17Af*. The total 
number of simulated haloes used for the analysis are 35381, 25418, 
7060, 5334, and 1840 in panels (a)-(e) respectively. 



analysis as the last subsection, and plot the results in Fig0| 
Since the LCDMc simulation has only one realization, no 
error bars are plotted for dP/dzf. However, because of the 
high mass resolution, the population of haloes within each 
mass range is so large that the random fluctuation among 
the bins is small. This also makes us be able to extend the 
analysis to massive haloes with up to 5000 member particles. 

Figure 0] continuously shows the change of the forma- 
tion redshift distribution with the halo mass. The pattern 
of the differences between simulation results and predictions 
is similar to that of small dark haloes. For the very massive 
haloes [FigUTe)] , the differences among the three models are 
small, and they agree with the simulation results reason- 
ably well. The results shown in Fig. |3f c) has mass range of 
0.18 — 0.36Af» which is almost the same as that of Fig.^Ja) 
(whose mass range is 0.17 — 0.35Af„), even though as used in 
these simulations are slightly different. Therefore, these two 
results can be compared for the same population of haloes 
with different mass resolutions. As can be seen, the results 
are consistent within the errorbars, though the population 
in the LCDMb simulation appears to form slightly later that 
that in the LCDMc. 



3.4 Resolution tests 

We consider halos with at least 100 particles, so the iden- 
tification of the halos is generally secure. There are still a 
couple of issues one should consider about the simulation 
resolutions. If the simulations are not started sufficiently 
early, some non-linear collapse could be missed or delayed, 
and the initial distribution of particles might still have an 
appreciable effect at the high redshift. Another worry is that 
some low mass haloes may be missed because of the force 
softening adopted in the simulations. To check these issues, 
we performed another two simulations that were run un- 
til z ~ 3.18. These two simulations have the same model 
parameters and the same initial fluctuations(including the 
phase) as the first realization of the LCDMb simulations, 
but one simulation was started at an earlier epoch z% — 72 
and the other adopts force softening r\ — 78 kpc that is twice 
large of the value used in the LCDMb simulation (39 kpc). 
The two simulations have also the data output at z ~ 4.78. 
In Figure 03 we plot the mass functions vf(v) of haloes at 
both redshifts for these two simulations, and compared them 
with that of the corresponding LCDMb realization, where 
v = [Sc(z) /a(m)} 2 . Except for very massive haloes that are 
rare objects (so their mass functions have large statistical 
fluctuations), the agreement of the mass functions among 
the three simulations is nearly perfect, especially for the 
mass range of 50 to 600 member particles that we are inter- 
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Figure 5. (a) The mass function of haloes at z ~ 3.18 in the three 
simulations: the dots for the simulation with Zj = 72, the triangles 
for the simulation with the force softening 78 kpc, and the open 
squares for the first LCDMb realization. The lines labeled with 
N p indicate the v values corresponding to 50 and 600 member 
particles respectively. The dashed line and solid line represent 
the predictions by the EPS model and the EC model, (b) The 
same as (a), but at z ~ 4.78. 



ested in this paper. Therefore we conclude that our results 
here are robust to the starting redshifts and to the force 
softening of the simulations. 



4 DISCUSSION AND SUMMARY 

This work was initially motivated by the fact that the 
semi-analytical models of galaxy formation which were con- 
structed based on the EPS theory for the merger tree predict 
systematically redder colours for dwarf galaxies than the ob- 
servations of faint blue galaxies (van den Bosch 2002b). It is 
known that the EPS theory fails accurately describing the 
mass function and the spatial correlation of small haloes 
(M <C M* , Jing 1998, Lee & Shandarin 1998, Sheth & Tor- 
men 1999). This difficulty motivated many theorists to re- 
formulate the formation of haloes, and reproduce the mass 
function as well as the spatial two-point correlation func- 
tion of haloes in numerical simulations (Sheth et al. 2001; 
Chiueh & Lee 2001, Lin et al. 2002). It is therefore interest- 
ing to test if the EPS theory, or the alternative models can 
correctly predict the distribution of the halo formation time 
(which is more closely related to the problem of faint blue 
galaxies) . 

With the help of a large set of high-resolution N-body 
simulations, we have found various degrees of success and 
failure for the analytical models at different halo mass. This 
is not unexpected since the excursion set approaches roughly 
model the growth of haloes but overlook some factors which 



could be important. For example, the models overlook the 
correlation between the fluctuations on different scales, do 
not include tidal striping, etc. If the conditional mass func- 
tion has problems, the merger trees constructed upon it may 
not correctly describe the growth of haloes (especially for 
small haloes) and therefore should be used with caution in 
galaxy formation models. We suggest that future studies of 
the conditional mass function should be pushed toward lower 
mass end using simulations with high mass resolution (like 
those used here). 

The implications of our results for the problem of the 
faint blue galaxies are quite clear. If the merger trees used by 
the semi-analytical models are replaced by the merger trees 
from the simulations, one would expect that the dwarf galax- 
ies will become generally redder, which contradicts with the 
observations more seriously. There are several ways to solve 
this discrepancy. First, the observed faint blue galaxies may 
not constitute a representative sample of small haloes. This 
could happen if many of red dwarf galaxies, because of their 
low star formation rate, have escaped from being detected. 
Second, star formation within small haloes is significantly 
delayed due to heating or re-ionization (see Mo& Mao 2002, 
and reference therein). Third, faint blue galaxies may expe- 
rience a recent interaction with nearby galaxies, which will 
trigger star formation. We have traced the trajectories of 
small haloes, and found that about 10% ~ 15% of the small 
haloes once pass the central part (with distance less than 
half of a virial radius) of a bigger halo which is at least 3 
times more massive than the small one. These haloes should 
have experienced a strong interaction. If gas in the small 
halo would not be stripped off, these strong interactions 
could trigger star formation and change the colour of the 
dwarf galaxies to blue. However these processes are so com- 
plex that more efforts are needed to work out the problem 
of dwarf galaxies. 

From the results of the LCDMb simulations, we found 
an earlier formation time for sub-M* haloes (about 10 12 so- 
lar mass) than the EPS prediction. This may indicate that 
galaxies are formed slightly earlier in the numerical simula- 
tions than in the EPS model. This result can explain a recent 
finding of Cimatti et al. (2002). Cimatti et al. measured the 
redshift distribution of galaxies with K < 20, and compared 
it with the predictions of semi-analytical models of galaxy 
formation based on the EPS merger trees. They found that 
the galaxies in the observation formed earlier than in the 
galaxy formation models of about 0.3 in redshift. Qualita- 
tively this can be seen in our Figure 3. Of course, the for- 
mation times of haloes and galaxies are not the same thing, 
but they should have a similar trend. 

There seems to be discrepancy between our results and 
those of van den Bosch (2002a) who found that the for- 
mation time distribution of the haloes in the GIF simu- 
lation deviates from the EPS prediction more severely at 
a larger mass. In fact, he checked for the haloes in two 
mass ranges: 5.6 x 10 11 < M < 1.1 x 10 12 fe _1 M and 
2.0 x 10 12 < M < 2.0 x lO 14 ^ 1 M Q (his Figure 4), cor- 
responding to the number of particles 40 < N p < 80 and 
140 < N p < 14000 respectively. The number of haloes in 
each mass bin is dominated by the lower mass end, es- 
pecially in the larger mass bin. For the small haloes of 
M < 1.1 x lO 12 /)." 1 M e , it is likely that his formation time is 
underestimated (i.e. formation delayed) because of the lim- 
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ited resolution, and so his agreement with the EPS theory 
becomes better. His result for large haloes should be com- 
pared with the middle panel of our Figure 3 (according to 
the effective mass). This comparison clearly shows that his 
results agree well with ours. 

In summary, none of the predictions by the EPS, EC 
and NCB model can match the simulation results at both 
low and high redshift. This discrepancy implies that the the- 
oretical mass conditional function needs further improve- 
ment. Our results may have an important impact on the 
models of galaxy formation. The "blue-color" problem of 
dwarf galaxies can hardly be solved solely by the gravita- 
tional interaction between dark matter. For sub-M, haloes, 
an increase of formation redshift by ~ 0.3 relative to the 
EPS model may relieve the difficulty of the semi-analytical 
models to explain the recent observation of K20 galaxies 
(Cimatti et al. 2002) at high redshift. 
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APPENDIX A: THE LINEAR GROWTH OF 
THE DENSITY PERTURBATION 



We assume that the universe is dominated by cold dark 
matter. The mass M of a halo is related with the co- 
moving length scale ro using a top-hat filter in real space 
M = (47r/3)7rroPo, where po is the co-moving mean mass 
density of the Universe. The variance of the density fluctu- 
ations, smoothed over a region of mass M, is given by 



a (ro) 



dk 
T 



k 3 P(k)W 2 (kr ) 



(Al) 



where W(kro) is the window function for top-hat filtering 
sin(fcro) cos(fcro 



W(kr ) = 3 , , 

(kro) 6 (kro) 2 

The power spectrum at present is given by 

P(k) oc kT 2 {k) 



(A2) 



(A3) 



which uses the self-similar primordial density spectrum with 
index n = 1. For the transfer function, we have used 

ln(l + 2.34g) 
v ' 2.34g 

[l + 3.89g+ (16.1g) 2 
and 
= k 

q ~ r^Mpc 



+ (5.46g) 3 + (6.7L?) 4 



-1/4 



(A4) 



(A5) 



(Bardccn ct al. 19861. In practice <r 2 (ro) was calculated up 
to an overall constant which is fixed by the choice of as = 
0-(8/i _1 Mpc). 

Collapsed halos are taken to be regions in the linear 
density field with the density contrast greater than some 
critical density contrast, S c . In practice we take into account 
of the linear growth by holding the variance a fixed and 
increasing the density thresholds at high redshift, 



* c (t(*)) = (l+*)»g£) * e . 

g{n(z)) 

The growth rate can be approximated as 



g(Q(z)) = -(I 
and 

O(z) = fi 



1 209^ fl 2 „ 

— i n h o 

70 140 140 



(l + 2 ) 3 



4/7 



i - fi m + (i + z) 3 n m 

for a flat A universe ^Carroll, Press fc Turner 1992) . 



(A6) 



(A7) 



(A8) 
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